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In this paper we first give a brief review of the variable cosmological constant model and its 
scalar field description. We mainly discuss two types of variable cosmological constant models: a 
power law and H power law models. A method to obtain all of the equivalent scalar field potentials 
and the effective equation of state of the two models is presented. In addition, the dynamics of 
such scalar field potentials and effective equation of state are discussed in detail. The parameters 
of the two models are constrained by current 307 high-quality "Union" SN la data set, baryon 
acoustic oscillation (BAO) measurement from the Sloan Digital Sky Survey (SDSS), 9 observational 
H(z) data derived from the Gemini Deep Deep Survey (GDDS) and the shift parameter of the 
cosmic microwave background (CMB) given by the three-year Wilkinson Microwave Anisotropy 
Probe (WMAP) observations. We also calculate and draw the picture of the Hubble parameter, the 
deceleration parameter and the matter density of the two models. Then, we show that the indices m 
and n in the two models have specific meaning in determining properties of the models. Moreover, 
The reasons that indices m and n may also influence the behavior of effective equation of state and 
scalar field potentials are presented. 
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I. INTRODUCTION 

In 1998, the discovery that the accelerated expansion of 
the Universe is driven by the dark energy (DE) from the 
type la supernovae (SN la) observations [lj greatly aston- 
ished the world. The Wilkinson Microwave Anisotropy 
Probe [2j, combined with more accurate SN la data [3| 
indicates that the Universe is almost spatially flat and 
the dark energy accounts for about 70% of the total con- 
tent of the Universe. However, we know little about the 
nature of dark energy except for its negative pressure. 
Therefore, a large number of works have been done in 
recent years to explain this mystery. 

The variable cosmological constant (hereafter VCC) [i| 
is one of the phenomcnological ways to explain the dark 
energy problem, because it is a straightforward modifi- 
cation of the cosmological constant A which enable itself 
to be compatible with observations. Looking back to the 
history, we can see that a lot of theorists have done nu- 
merous works to search for the theoretical foundation of 
the VCC models and also investigate the properties of 
the VCC models @. In [f|, a model of A cx a~ 2 was pro- 
posed, requiring that the cosmic density p would equal to 
the Einstein-de Sitter critical density p c , which leads to 
a closed Universe, without singularity, horizon, entropy 
and monopoly problems @. After that, it was also sug- 
gested a model A cx a~ 2 (A should be independent of 
fi) with different initial conditions by Q, which firstly 
pointed out that time-dependent A leads to the creation 
of matter or radiation. Besides, a lot of work were done 
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to propose straightforward models relating A to the Hub- 
ble parameter H(z): A cx H 2 @, S 0, ■ Furthermore, 
people also constructed a large number of phenomenal 
VCC models to describe the dynamics of the Universe 
and there is a list in Ref. [i| summarizing the proposed 
models [l2j . There are also several papers concerning the 
observational constraints about the VCC models 

In addition to the VCC models, scalar fields such 
as "quintessence" [l4j], "phantom" jlj| and "quintom" 
have been introduced to effectively describe the dy- 
namic dark energy, which are distinguished by the ef- 
fective equation of state (hereafter EEoS): wde > — 1, 
Wde < — 1 and wde across —1 respectively. These mod- 
els are inspired by the fact that a decaying vacuum en- 
ergy which has the very high energy density at early time 
should be sufficiently small at present to meet the current 
observation requirement, so they should evolve dynami- 
cally. In order to obtain the corresponding quintessence 
potentials, the reconstruction equations were derived and 
addressed the feasibility of the approach by Monte-Carlo 
simulation [l^] ; it was also constructed the general scalar- 
field dark energy model and developed a method 
to construct them directly from EEoS function w<j,{z) 
[l^, |2(j. Moreover, some works have been done to recon- 
struct the scalar potential from the scalar-tensor theory 
and investigate the modified Newton theory [2l| . 

As a major part of our work, we analyze the EEoS 
and reconstruct the potentials for two main types of 

VCC models the a power law and the H power law 

models from the point of view of dynamic scalar fields. 

This work is necessary for people who are interested in 
the coupled dark energy and dark matter [22| . because 
such models may avoid a lot of realistic problems such 
as the coincidence problem. In addition, it is discussed 
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how such phenomenological models can be explained as 
a classical scalar field decaying into a perfect fluid which 
might be interested by those who want to search for the 
gravitational theory other than the general relativity, be- 
cause the Lagrangian in VCC should be different from the 
Einstein-Hibert action in general relativity. Thus, this 
part should be essential for people to see the possible 
forms of VCC models and its corresponding scalar fields 
which are expected from string theory or supergravity. 

Another main part of this paper is to give an obser- 
vational constraint on the VCC models and explain the 
properties of cosmological parameters. This part of work 
is the basic analysis to determine the right form of VCC 
models from the observational requirement. 

This paper is organized as follows: in section 2, we 
search for the EEoS, the reconstruction equation and ef- 
fective potentials for all the a power law models, general- 
izing previous work from Ref. [231 ]. Next, we use the 
current observational data, including 307 high-quality 
"Union" SN la data set, baryon acoustic oscillation from 
SDSS, 9 observational H(z) data and CMB shift parame- 
ter from WMAP three years result, to constrain the index 
in this model. In addition, we analyze the properties of 
the dark energy density, the dark matter density and the 
deceleration parameter in this model. In section 3, par- 
allelling to section 2, we generalize the work from Ref. 
[HI and discuss the EEoS, the reconstruction equation 
and reconstructed potentials for all of the H power law 
models. In addition, we also give one example of this 
type of model to prove the effectiveness of our method. 
Then, we use data to constrain the cosmological param- 
eters of this type of models and analyze the properties 
of the dark matter density, Hubble parameter and the 
deceleration parameter. The concluding remark will be 
presented in the last section. 

II. a POWER LAW MODELS AND 
CORRESPONDING POTENTIALS 

For a generalized VCC related to the scale factor a, we 
can write 

A = Ba~ m , (1) 

where B is a constant with the dimension of mass square 
and we call it the dark energy amplitude. We assume 
the VCC is proportional to the scale factor power — m, 
and power index m plays a significant role in determining 
the dark energy behavior as discussed below. Then, the 
dark energy density and the Friedmann equation can be 
written as (Note that a = 1, ^ = -H(l + z)-^.) 

p<£\z)= AM^ = B(l + zrM^. (2) 

3A%H a = p$+p < £>. (3) 

For simplicity of calculations we assume spatial flatness 
(k = 0) which is motivated by theoretical considerations, 



such as inflation, and also confirmed by current obser- 
vations such as WMAP three years result Q. Our re- 
sults can be easily generalized to the case with a spatial 
curvature. We denote M p i = (8ttG)~^ as the reduced 
Planck mass and use superscript (a) here to denote a 
power law model, we will also use superscript (H) to de- 
note H power law model in the next section, pm is the 
dust matter density with the present value 

pH = 3F 2 M p 2 ,O ro0 . (4) 

Thus, we get 

£ = 3# 2 (1-Q m0 ). (5) 

In our practice, there is clearly only one degree of freedom 
in the a power law model, which is the power index m. 
As for the VCC models, it is rather natural to consider 
the interaction between the dark matter and dark energy 
[24l ]. Therefore, we should introduce an interacting term 
Q(z) with 

pV>+3Hp$=Q{z), (6) 

p^+3H(p^+p^) = ~Q(z), (7) 

and the total energy conservation equation 

ptot + 3H(p tot + ptot) = 0, (8) 

still holds. Since VCC is the generalized from of cosmo- 
logical constant, so it satisfies p^ = —p^- The Eq. ([7]) 
leads to 

Q(z) = -p^ = A(l + z) m H, (9) 

where 

A = 3F 2 (1 - n m0 )mM^, (10) 

which means that the interaction is explicitly determined 
only by the evolution of dark energy density. 

A. Interacting dark energy and reconstructed 
potentials in the a power law models 

The anterior Eqs. ([6]) and ([7]) are the standard in- 
teracting dark energy equations and the function Q(z) 
represents the interaction between dark energy and dark 
matter. Since the interaction may not be directly observ- 
able, it is interesting to search for the phenomenologically 
equivalent potentials which encode some properties of the 
interaction. 

One way to search for such a theory is to express 
the VCC models in a field theory language, so the most 
straightforward way might be the scalar field description. 
If one could find such a description of the VCC models, it 
is natural to extend the scalar field description to other 
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space-time and other gravitational theory like superstring 
theory. In addition, this description is very useful since 
it provides a path to quantize the scalar field, which can 
help people to understand the fundamental theory of the 
phenomenological VCC models. Furthermore, the proce- 
dure to obtain a scalar field description of a phenomeno- 
logical model could be applied to other models. 

From this point of view, we want to see what are the 
EEoS and dark energy potentials in the a power law mod- 
els. Changing the form of Eq. ([7]), we have 

p^+3H(p^+p^ + ^)=0, (11) 

so we could see the interaction Q(z) contribute to the 
effective pressure 



(a) _ (a) Q(z) 

Peff - P\ + 3H 

(a) , Q(z) 



= -Pa 

so the EEoS of dark energy is 



.(<*) 
,J eff 



P ( a\z) 



3H 



to 



(12) 



(13) 



so we obtain this result from the point of view of in- 
teracting dark energy [2(| HH . In Eq. ([IB]) , the power 
index to is a constant, so the EEoS in the a power law 
models are all constants. We will see in the following sub- 
section that the best-fit of index to constrained by cur- 
rent combined observational data is —0.09, which means 
< —1, so VCC is phantom-like [HI]. However, for 2a 



,(«) 

J eff 



confidence level we cannot rule out the possibility that 
to > (quintessence-like), so we should consider both 
the phantom and quintessence scalar field potentials for 
the VCC model. 

For a spatially homogeneous and isotropic scalar field, 

the effective energy density p^ and pressure V^eff can ^ e 
written as 



i . 



(14) 



(15) 



phantom (quintessence) scalar field, and V^j((f>) is the 
effective scalar field potential for the a power law mod- 
els. At the same time, the effective energy density pj^ 
and pressure P^tf are given by the interacting dark en- 
ergy equations 



(a) 

Pa 



(a) 

-Pa 



Q[z) 



3H 



A 2 (i + Z y 



(16) 
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where 



and 



rn s 



A 2 = A 1 {-1 + -). 



(18) 
(19) 



We define the dimensionless quantities 

I = <f>/M pl , V eff = V eff /ZHlM 2 pl . (20) 

Thus, the scalar field potential can be written as a func- 
tion of redshift z 



where 



V {a) (7\ - -(n {a) - n {a) ) 
= A,(l + z) m . 



A 3 = (i-n m0 )(i--). 
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Combining Eqs. IT4]) and (fTS")) . we have 

C x 



-r = T 



where 



dz (1 + Z )[C 2 + C 3 (l + z) 3 ' 



d = (l-r» m0 ) 5 x [H(3-to)]" 



C 2 = 3(1 - O m0 ), C 3 = {3Q m0 - to). 



(21) 
(22) 

(23) 

(24) 
(25) 



The upper (lower) sign in Eq. |23|) represent <j> > 
(4> < 0). In fact, the sigh is arbitrarily determined by 
assumption, as it can be changed by (f> — > —(f). We choose 
the upper sign in the following discussion. If we shift 
(j>o value, the potential in the following figure will be 
shifted horizontally, but the shift doesn't influence the 
whole shape of the potential. The field could be inte- 
grated analytically as 



(z) = Ca x tanh 1 



respectively, where upper (lower) sign represents the where 



91 i 

(C 2 + C 3 (l + z) 3 - m )i 1 



Cx = 



V3 L 3 



(26) 



(27) 



We let the integral constant equals to zero since the initial 
value of field is meaningless. Solving this for (1 + z) and 
substituting the result into (|2Tj) , we obtain the potential 
of the a power law model 



V e ( ;)W = A 3 [^ x (coth a (^) - 1)]**. (28) 
We use the best fits value for to in the next section to 
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FIG. 1: Reconstructed potentials for a power law model. Here 
we set fl m o — 0.28. 



draw the pictures for the a power law models' equiva- 
lent potentials. There are three main characteristics for 
these potentials: First, they are all runaway type and 
the whole shape doesn't change if <fi is shifted horizon- 
tally. Second, Eqs. (f2Tj) and (126)) determine that cj) in- 
creases and V^j(4>) increases as the redshift z decreases 
from large value to —1, which means that the dark en- 
ergy potentials increase as the Universe expands. From 
the figure, we could see that more negative the value of 
to is, the sharper the potentials will increase as the field 
evolves. One of the possible explanation for this phe- 
nomenon is that the more negative value of m means the 
dark matter "decays" into dark energy quicker, so the 
dark energy increase its potential value and energy den- 
sity faster. Third, the hyperbolic coth function in the 
expression (f2"5f makes the a power law potentials have 
the asymptotic value. This is very interesting because 
one could obtain such behavior in general in the super- 
symmetric QFT. This runaway form of potential is also 
the one expected in the unstable D-brane system in su- 
perstring theory (26j. 



B. Hubble parameter and results of the constraints 



In this subsection, we want to constrain the parameter 
to from combined observational data, so we should obtain 
the Hubble parameter and the luminosity distance. We 
can change the variable t to redshift z in the Eq. © to 
figure out the analytical expression for the matter density 

p%\z; to) = 3HM [^(1 + zf + D 2 (l + z) m ], (29) 



where 



C 3 



-, D 2 



3 — to 



(1 — fl m o). 



(30) 



We can easily note that the effect of VCC is just like a 
small perturbation to the evolution of the matter density. 
If the evolution behavior of the matter density doesn't 
deviate much from (1 + z) 3 behavior (ACDM), the value 
of to should be very near zero, which indicate that even 
though the dark energy is not constant through the evolu- 
tion of the Universe, it at least should evolve very slow. 
This property will be convinced through observational 
constraints on power index to in the following subsec- 
tion. This equation is essential for our purpose to solve 
the Hubble parameter in the following subsection. 

As there is only one free parameter in this kind of 
power law models, it is rather easy to obtain the best fit 
value from the current observational data. We do this 
fitting using the high quality type la supernovae, baryon 
acoustic oscillation from SDSS, observational H (z) data 
and the the shift parameter of the cosmic microwave 
background (CMB) given by WMAP three years results. 
We discuss this problem in the framework of interacting 
dark energy and accelerating Universe (see relevant work 

Integrating the Eq. (J3J), we have the following equation 



H(z) = ff [£>i(l + z? + D 3 (l + z) m ]~^ (31) 



where 



D,= 



C-2 



(32) 



Then we integrate the Eq. (|31j) and follow the guideline 
in Appendix (B) to obtain the x 2 formula to do numeri- 
cal fitting. Our result is much tighter than the previous 
fitting results [27| due to the more precise data we use. 
Under the combined data sets SN+BAO+OHD+CMB 
constraints, the 3a values for the power index to are 



777 = -0.09 



+0.08+0.12+0.19 
-0.11-0.20-0.29- 



(33) 



From the above results, we are still not able to rule out 
the positive to value within 2er confidence level, so we 
need more precise data to constrain the VCC models' 
power index in the future. However, no matter whether 
the to is negative or positive, it is very near zero, suggest- 
ing that even though the dark energy is not a constant, 
it should evolve very slow. The results for the best-fit 
and lcr values are shown in Table 1. From our results, 
the f2 m o is always around 0.28, which is consistent with 
the WMAP three years results . 



Matter density and deceleration parameter of 
the a power law models 



Having the matter density Eq. (|2"9"|) and the confidence 
region of parameter to, we can plot the matter density 
and dark energy density as a function of redshift z as 
FIG. [3]and FIG. HJshows. We put the curve representing 
the standard matter density equation in ACDM model 
for comparison. 
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Models 


SN+BAO 


SN+BAO+OHD 


SN±BAO±OHD±CMB 


a power 
law 


m 


-o.i2t° :« 


-0 1Q+ - 29 


-o.o9±°:™ 




0.28i° ;° ! 


n QQ+O. 04 


29+ 03 


H power 
law 


n 


-n 26+ - 67 

u - ZD -0.74 




-0 TS+ 014 




0.28 ±0.04 


0.27 ±0.03 


0.29 ± 0.03 



TABLE I: Results of the fitting for the two models. 
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FIG. 2: Contour map for the parameter m versus f2 m o in the 
a power law model. Green dashed lines represent SN±BAO, 
blue lines represent SN±BAO±OHD, and red lines represent 
SN±BAO±OHD±CMB. 



From FIG. [3] and dj if to is positive, the dark energy 
density changes and slower the matter density decreases 
as a function of redshift z. That means the dark energy 
" decays" into dark matter field and makes the matter di- 
lute more slowly with the cosmic expansion, vice versa. 
As a result, the dark energy density will continuously de- 
crease as the Universe evolves. On the contrary, if to < 0, 
the dark matter transforms into dark energy and the dark 
matter energy density decreases more sharply than the 
usual (1 ± z) 3 behavior (ACDM). Then the dark energy 
increases its energy density and realizes the "Big Rip" 
in the future due to this matter changes, so effectively it 
resembles the phantom dark energy (see Eq. ([TBf ). 

It is easy to see this property of index m through the 
"decay rate" e in Ref. [29j . The "decay rate" e is defined 
as the matter density's deviation from the standard evo- 
lution, i.e., 

Pm = p m0 a _3+£ , (34) 
where p m o is the current matter density. In our case, e is 




FIG. 3: Matter density for a power law models. Here we set 
£2m0 = 0.28. 
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FIG. 4: Dark energy density for a power law models. Here 
we set Qmo = 0.28. 



not a constant but a function of redshift z. In addition, it 
is straightforward to verify that m is the index to distin- 
guish the sign of e, as to > 0, e(z) > 0; vice versa. Thus, 
the relationship between index to and the "decay rate" 
e represents whether the dark energy " decays" into dark 
matter or the inverse. Moreover, it is also easy to con- 
firm that the value of e(z) is generally compatible with 
the confidence region provided by [29| . 

Having obtained some meaning of the index m and 
its confidence region in the a power law models, we can 
directly find the evolution of the deceleration parameter 
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in this kind of model. 



where 



(36) 



From FIG. [5j we can understand the following charac- 
0.5 F I - 



N 




results are rather consistent with this work because this 
best fits value for zt will lead to the best fits region for m 
[—0.43, 0.16], and our 3tr results for m is just within this 
region, suggesting that our constraints on m and £l m o are 
very tight. At the same time, our results are also very 
consistent with other relevant constraints on the " equiv- 
alent" redshift when p m {z eq ) = p\(z eq ) [3l| . 

Therefore, in this subsection we conclude that the 
power index m of the a power law model is not only asso- 
ciated with the dark energy density, but also a meaning- 
ful index to determine whether the dark matter " decays" 
into dark energy or the inverse. Moreover, it determines 
the "decay rate" e, i.e., the intensity by which dark mat- 
ter changes into dark energy. Meanwhile, it affects the 
deceleration parameter of the Universe and the transi- 
tion redshift zt when the Universe was changing from 
deceleration to acceleration. 



III. H POWER LAW MODEL AND ITS 
RECONSTRUCTED POTENTIALS 
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FIG. 5: Deceleration parameter for a power law models. Here 
we set Qmo = 0.28. 

tcristics about the deceleration parameter in the a power 
law models: First, different a power law behaviors (in- 
cluding ACDM) have common deceleration parameters 
at present time q(0) = |fi m o — 1- Second, they have dif- 
ferent values of redshift z when the Universe was chang- 
ing from deceleration to acceleration, so the transition 
redshift zt for different models are determined by the 
following equation 



ZT = [ 



;2-m)C 2l _L. 



C 3 



1. 



(37) 



which can be determined by the results of the con- 
straints. In addition, we could calculate the different 
zt corresponding to m as the best fit, ler, 2<r confi- 
dence values and the ACDM: z T (m = -0.29) = 0.57, 
z T {m = -0.20) = 0.61, z T (m = -0.09) = 0.67, z T (m = 
-0.01) = 0.72, z T (m = 0.03) = 0.75. So the larger m 
is, the earlier the Universe changes from deceleration to 
acceleration. Thus, we obtain this transition redshift zt 
from the general a power law form and the results should 
be applicable for all of the specific power law behaviors 
[2{|. Third, the more negative value of m is, i.e., more 
sharply dark energy density changes, the faster Universe 
accelerates. This result is also compatible with the FIG. 
3 in [H. 

Furthermore, we can also use the constraints on 
the deceleration parameter and transition redshift to 
see whether our results are consistent with relevant 
constraints [30L l3ll . |32| . In Ref. [30], it shows that the 
best fits for the transition redshift is zt — 0.781q27. O ur 



In this section, we will discuss another type of VCC 

models A is associated with Hubble parameter H — 

which is an important type presented in Ref. 0]. 

In this type of model, the VCC can be written as 



A = CH V 



(38) 



where C is a constant with the dimension of mass 2 — n, n 
is the only parameter in this kind of models which needs 
to be fitted by observational data. Then, the dark energy 
density and the Friedmann equation in this model can be 
given by 



Pa 



AM; 



CH n M 2 ph 



3M 2 ff 2 = pCH) 



(39) 



(40) 



The amplitude C is determined by the current value of 
matter density and the Hubble constant 



C* = 3i^-"(l - fi, 



(41) 



Then, we consider that the VCC indicates that there is 
an interaction between the dark matter and dark energy. 
Therefore, let's assume that dark energy and matter ex- 
change pressure through the interaction term W(z) with 



P, 



H) 



3Hp 



W(z), 



.(H) 

Pa 



mpT ] 



-W(z) 



(42) 



(43) 



which maintains the total energy conservation equation 
ptot + SH{p to t +Ptot) — 0. Since the VCC is the general- 
ization of the cosmological constant, so it should satisfy 
then Eq. |g3J) leads to 



(H) (H 

Pa = -Pa 



W{z) 



.(H) 

-Pa 

G 1 H n {l + z)H'(z), 



(44) 
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where 



G x =3nH^ n {l~n m0 )M^. 



(45) 



for dark energy density and pressure through definition 
([55)) and the interacting dark energy Eq. ([4"3")l . 



Pa 



G, 



(52) 



A. Interacting dark energy and reconstructed 
potentials in the H power law model 

In this subsection, we want to see the potential that 
dark energy mimics the VCC. Although the interaction 
between dark energy and dark matter might not be di- 
rectly observable, the effective potential could encode 
some information about the interaction. Thus, we are 
looking forward to solving the EEoS and reconstruct the 
dark energy potentials of the H power law models from 
the standard interacting dark energy Eqs. (|42|) and (|43|) . 

Transforming Eq. (|4"5|) , we have 

pP+W(pF+pP + ^) = 0, (46) 

from which the interaction changes the effective pressure 
of this model, i.e. 



AH) _ IB) W(z) 
Peff - Pa + 3H > 



so the EEoS of dark energy is 



, Peff n G 2 (z) 

^ a ff - — — 7777 — — J- 



ft p W ' ' 2G 2 (z) + G 3 ' 



(47) 



(48) 



where 



G 2 (z) = n m0 (l + z) 3 ^-i n \G 3 = (1 - Slmo). (49) 

These EEoS are functions of redshift z, in contrast to 
the a power law models, where the EEoS is constant in 
the Eq. (flU]) . It is quite interesting that the sign of in- 
dex n also determines whether this dark energy likes the 
quintessence or phantom. Moreover, this type of EEoS 
is affected by the value of fl m o, while the EEoS in the a 
power law models are not. Since the constraint from the 
current observational data suggests that the best-fit for n 
is negative but cannot rule out the possibility of positive 
constant n (see discussion in next section), we construct 
the potentials for the two cases. The energy density and 
pressure density of the quintessence field for this model 



T^ 2 + K ( /?(0) 



T 2< 



(H) 

Pa , 



v {H) 

Peff 



(50) 



(51) 



where the upper (lower) sign represents the phantom 
(quintessence) dark energy, corresponding to n < 
(n > 0). At the same time, we can obtain the expressions 



(H) 
Pa 



- -A 1 [G 2 {z) + G d f^ 1 
x [(l-£)G 2 (z) + G 8 ]. 



(53) 



Then, the effective scalar potential can be written as a 
function of redshift z 



V e f/(z) = G 3 [G 2 (z) + G 3 ]^ 



x [ ( l--)G 2 (z) + G 3 ] 



(54) 



Using the Eqs. (|5D|) and (p)Tj) . we can obtain the differ- 
ential form of scalar field 



where 



# _ G 4 (l + z)*-*' 
~d~z~ T G 2 (z) + G 3 



G 4 = (-H)i(O ro0 (l-O m0 ))3 



(55) 



(56) 



In general, Eq (|55|) could be solved analytically so we 
obtain the following field equation 



(z) = 0o T G 5 arctan[G 6 (l + z)^ 1 '^], (57) 



where 



G 5 



2 ,2 i Qmo a /rnX 
"(o l"l) 2 :G 6 = (- p: ) 2 j (58) 



2 - n 3 



1 — SimO 



and the upper(lower) sign applies if > (0 < 0). In 
fact, the sigh is arbitrarily determined by assumption, as 
it can be changed by 4> — > — cj>. We substitute (1 + z) for 
<t> into Eq. (|54"|) to get the result of potential V^f)(4>). 



v e f f ) (4>) = G!- n [i + t a xi 2 C- 



G 



(59) 



We give the following three examples of phantom poten- 
tials for this kind of models (see FIG. [6]). As is shown in 
FIG. [HI the effective phantom potentials also have some 
characteristics: For one thing, they are all runaway type 
potentials and even if we change the initial value <fo, 
the curves shift horizontally with the whole shape un- 
changed. For another thing, the meanings of the poten- 
tials are clear: as the Universe is expanding, the value 
of 4> becomes large and the field slowly rolls upon the 
potential, which makes the EEoS very close to —1. At 
the same time, the dark matter field gradually "decays" 
into dark energy, so the dark energy density increases its 
energy density as the Universe expands. 
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0.0 
0.00 



0.20 



FIG. 6: Reconstructed potentials for H power law model. 
Here we set ttmo = 0.28. 



B. One specific examples of reconstructed 
potentials for the H power law models 

In Ref. [H, there is a list of proposed H power law 
models proposed by different authors through various 
perspectives. In addition, Ref. [23[ gives two examples of 
scalar potentials from the point of view of the scalar field 
description. In order to show the effectiveness of our re- 
construction, we derive one analytic results of the scalar 
potentials using the methods in previous subsection. 

n = 2 is an interesting Q case of the VCC model 
and investigated by many authors 0, [23| . We substitute 
n = 2 into Eq. (fB"4")) to obtain 

' " (60) 



V${z) = (1 - Q m0 )(l ~ |^mo)(l + z) 
and the field ([55]) could be integrated as 
0(a) = o -3G 7 hi(l + *), 
where </>o is the initial value of field <fi and 
G~i = f2 m o(l — ^mo)- 



3Q„ 



(61) 



(62) 



Thus, we obtain the potential by substituting (1 + z) for 
field A 



"eff 

where a 



KV7W = ( 1 -^o)(i 



i 



- Q r 



-ot(4>-cj>o) 



(63) 



30„ 



-) 2 . This form is rather consistent with 



that in Ref. [23| , which demonstrates the effectiveness of 
the reconstructing method in this paper. This potential 
is one of the simplest runaway types which represents 
the particle creation in the phenomenological VCC mod- 
els so it could be interpreted as some kind of "coupled 
quintessence" [22j ]. Meanwhile, it is also easy to see that 
all the VCC potentials are associated with the exponen- 
tial function, which leads to its runaway behavior, indi- 
cating that they might be easily obtained in supergravity 
and unstable D-brane systems [26| . 



C. Hubble parameter and results of the constraints 

on n 



From Eqs. (gOJ), S3) and (g3]), we can obtain the dif- 
ferential equation for the Hubble parameter 



H (*)- 



1 



-(3H - K^' 1 ) =0, 



2(1 + *) 
where 

The Eq. ([64]) can be solved analytically 
H(z)=H Q [G 2 (z) + G 3 ]^ 



(64) 



(65) 



(66) 



Then, we follow the procedure in Appendix (A) and 
(B) to do the numerical fitting and finally we can find 
la results as Table 1 shows. Under the combined 
SN+BAO+OHD+CMB constraints, the best fit, la, 2a 
and 3a values for parameter n is 



„ _ n 1 c+0. 14+0. 23+0. 25 
a — U.lO_ 17 _ .26-0.43- 



(67) 



From the result (|67[) . no matter whether n is negative 



-1 - 



-2 - 









SNIa+BAO 








1 SNIa+BAO+OHD+CMB 




/ SNIa+BAO+OHD 









0.0 



0.2 



0.4 



0.6 



0.8 



1.0 



a, 



FIG. 7: Contour map for the parameter n versus f2 m o in 
the H power law model. Green lines represent SN+BAO, 
blue lines represent SN+BAO+OHD, and red lines represent 
SN+BAO+OHD+CMB. 

or positive, it always very near 0, indicating the slow 
evolution of dark energy. We could plot the EEoS PS]) 
as a function of redshift z and compare them with other 
models [28| and observational results [33[ . This EEoS has 
three major properties: For one thing, the confidence re- 
gion of this type of models is mildly consistent with the 
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FIG. 8: The redshift dependence of the EEoS for H power 
law models. Here we set fl m o = 0.28. 



FIG. 9: Matter density for H power law models. Here we set 
£2 m o = 0.28. 



results which were obtained by using CMB and Clus- 
ters data [33| . indicating that it is a competitive model 
waiting for the examination by future observations. For 
another thing, within la confidence region, the EEoS is 
slightly less than —1, which implies that it resembles the 
phantom field. In contrast, within 2a region it is possible 
that the EEoS is greater than —1, so we cannot rule out 
the possibility that the dark energy is quintessence like. 
Thus, this type of models could really represent a large 
kind of dark energy models phenomenologically. Further 
more, if z become larger, all of the EEoS in these mod- 
els have their own asymptotically constant value, which 
is rather similar to that of the "quiessence model" [34| . 
The constant EEoS means that the proportion of kinetic 
energy to potential energy is constant. Thus, the whole 
dark energy density increases or decreases, suggesting 
that the VCC models corresponds to a dissipative sys- 
tem of dark energy [f| 0, • When redshift z approaches 
— 1, all of the EEoS w e ff(z) approach —1, indicating the 
Universe will enter the de-Sitter phase in the future. 



n < 0) decreases more sharply than the standard (1 + z) 3 
behavior (ACDM). As a result, the parameter n is not 
only a power index of the H power law models, but also 
an important signature to distinguish whether dark en- 
ergy "decays" into dark matter or the inverse process, 
just as the index m in the a power law case. Further 
more, comparing FIG. O and FIG. [21 we can discover 
that the two graphs are very similar to each other, which 

means that the two types of VCC models the a power 

law and the H power law models really share some 

common features if the parameters are all constrained by 
observational data. This could be understood as follows: 
both the scale factor a(t) and the Hubble parameter H{z) 
describe the evolution of the Universe; if the Universe is 
expanding canonically, a(t) will increase while the H(z) 
will decrease, so the difference between the two types 
might only lie in the sign of the power index. 

However, if we plot the Hubble parameter by selecting 
some ideal value of n, we could see that they indicate the 
different fates of the Universe. We see from Fig. (jTUJ) 



D. Matter density, Hubble parameter and 
deceleration parameter of the H power law models 

From the Hubble parameter Eq. (|64[) . we can obtain 
the matter density 

P W(z\n) = mlM 2 pl G 2 (z)[G 2 {z) + G 3 }^ . (68) 

Then, we plot the dark matter density and Hubble pa- 
rameter of the H power law model. In order to compare 
with ACDM model, we also plot n = curve in one 
graph. Thus, FIG. [9] helps us to analyze the properties 
of the matter density in this model: For one thing, we can 
see that for the positive n, the dark energy density de- 
creases and the matter density dilutes more slowly as the 
Universe evolves, vice versa. On the contrary, n < rep- 
resents the dark matter changes into dark energy because 
the dark energy density increases as time evolves and the 
corresponding curve (for example, the curve in the graph 
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FIG. 10: Hubble parameter for H power law models. Here 
we set Qmo = 0.28 and H(0) = 72km/s/Mpc. 

that for any selected value of n, the Hubble parameter 
cannot deverge, so the H power law model does not in- 
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dicate the "Big Rip" phase in the future. This property 
could be understood from Fig. ([5]), because whatever the 
current value of w e ff is, they will all tend to be — 1, so 
the Universe will enter a de-Sitter phase in the future. 
This property is rather different from a power law model 
since u> e // is a constant in that model, so the dark energy 
will be always a phantom if m < 0, which surely results 
in a phantom Universe with the "Big Rip" phase in the 
future (see Fig. (g])). 

Having obtained the matter density, we can derive the 
deceleration parameter in this model 



,(H) 



1 2 - K 2 (z) 
' 21 + K 2 (z)' 



where 



K 2 (z) 



9. 



mO 



1 — ^mO 



(1+*) 



3(1- 



(69) 



(70) 



Then, we can plot this deceleration parameter and also 
compare the curves with that of ACDM. From FIG. ([IT]). 



N 




FIG. 11: Deceleration parameter for H power law models. 
Here we set Q. m o = 0.28. 

we could analyze the properties of the deceleration pa- 
rameter in this H power law model. For one thing, just 
as the a power law case, all of the H power law models 
share the same deceleration parameter at present time, 
which is q(0) — §17 TO o — 1, but they also share the same 
deceleration parameter values in the future (a — > oo) 
as q(z = —1) = — 1, which is different from that of 
the a power law model. For another thing, transition 
redshift zt varies differently according to the different 
curves, which indicates that different values of n affect 
the expansion of the Universe distinctively. The transi- 
tion redshift zt for different models are determined by 
the following equation 



ZT = [ 2(i-n mo)] ^_ h 

"mO 



(71) 



We could calculate its value with respect to different in- 
dex n numerically (see Table 2). It is also easy to see 
that the value of zt within la confidence region is gener- 
ally compatible with the result in Ref. (29|. In addition, 



from Fig. ([TT|) . the transition redshift is in the range 
[0.57,0.77], which is just in the best fits region of tran- 
sition redshift in Ref. [3(| and [3l|, indicating that our 
numerical constraints on the parameters n has been al- 
ready very tight compared with relevant work [30l l3ll. 132]] . 
Moreover, different values of n correspond to different 
curves with distinctive shapes. For one thing, within lcr 
region, n is definitely negative (such as n = —0.15), this 
makes the dark energy density increases and dark mat- 
ter "decays" into dark energy as time evolves. The more 
negative the index n is, the more quickly the dark energy 
density p A increases and the faster the Universe accel- 
erates. For another thing, within 2a confidence region, 
there is a certain probability that the index n is positive, 
which indicates that the density of dark energy is de- 
creasing, so the acceleration is relatively small compared 
to the negative value of n. To conclude this subsection, 
we derive the evolution of the matter density, Hubble pa- 
rameter and the deceleration parameter in the H power 
law models. We note that the index n is just like the 
index m in the a power law models, which not only re- 
flects whether the dark energy "decays" into dark matter, 

but also affects the acceleration of the Universe the 

deceleration parameter q(z) and transition redshift zt- 



IV. CONCLUDING REMARKS 

In this paper, we develop a method to reconstruct po- 
tentials in the VCC models directly from the definition 
of the energy density and pressure of the scalar field. 
We also give one example of the reconstruction for the 
H power law models. First, these potentials have some 
relationship with the exponential function, which is ex- 
pected in supersymmetry theory and unstable D-brane 
system in superstring theory [26j. Second, as the Uni- 
verse is expanding, the value of cj) becomes large and the 
field slowly rolls upon the potential. At the same time, 
the dark matter field gradually "decays" into dark en- 
ergy, so the matter density dilutes more sharply than the 
standard (1 + z) 3 behavior. It is worth noticing that 
the reconstruction equations presented here are not lim- 
ited to searching for the scalar field description of such 
phenomenological VCC models. Generally, it could give 
people the possibility to find the scalar field versions of 
other phenomenological models and even quantum grav- 
ity models such as holographic models [35j and vacuum 
fluctuation model [3t| . 

We also investigate constraints on the VCC models — 
- the a power law and the H power law models — 
from current combined cosmological observations, using 
high quality supernovae, baryon acoustic oscillation from 
SDSS, observational H(z) data derived from Gemini Deep 
Deep Survey (GDDS) and the CMB shift parameter from 
WMAP three years result. We consider a spatially flat 
FRW Universe with matter component and VCC compo- 
nent. For the VCC models, such as the a power law and 
the H power law models, the power indices m and n play 
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J 



n 


-0.41(2cr) 


-0.32(1<t) 


-0.15(the best fit) 


-0.01(la) 


O(ACDM) 0.08(2cr) 




0.57 


0.60 


0.66 


0.72 


0.73 0.77 



TABLE II: transition redshift zt in H power law model. 



a very significant role in determining the evolutionary 

behavior of the space-time as well as the ultimate fate of 

the Universe. According to the combined constraints, the 

best fit, Id, 2(7 and 3c values for the indices to and n are 
_ n no+o.08+0.12+0.19 _ _ n -, r+o. 14+0.23+0. 25 



m 



re- 



0.11-0.20-0.29' " ~~ "-""-"-0. 17-0. 26-0. 43 

spectively. These results cannot fix the value of to andn, 
but at least indicate even though the dark energy may 
change with time, it will evolve very slow since the value 
of to andn are always very near zero. The indices of VCC 
models suggest the interaction between dark energy and 
dark matter: First, to > and n > represent that 
the dark energy "decays" into dark matter, while m < 
and n < represent the inverse process. The more neg- 
ative the indices to and n are, the faster such transitions 
happen. Second, the indices to and n affect the decel- 
eration parameter and the transition redshift Zt- The 
more negative the indices to and n are, the more portion 
dark energy takes in the whole Universe budget in the 
future due to the matter "decays"; thus, faster the Uni- 
verse will be accelerating. Also, the more negative the 
indices to and n are, the less portion they took into the 
whole Universe budget in the past and the smaller value 
of transition redshift zt is. Therefore, the indices m and 
n are the important signatures to judge whether the Uni- 
verse accelerates more drastically than the ACDM model. 
Third, the indices m and n are also the essential indica- 
tors to understand the properties of VCC. For one thing, 
the best fit values for to and n suggest that the EEoS 
of dark energy are real numbers or functions of redshift 
z at the region [—1.5,-1], indicating that the dynamic 
scalar fields of VCC are phantom-like. For another thing, 
there are still some probabilities for the a power law and 
the H power law models that the dynamic scalar field 
of VCC is quintessence-like, so the VCC models are the 
phenomenal models representing a variety of other dy- 
namic dark energy models. Moreover, the EEoS of the 
a power law models are constant, while those of the H 
power law models are functions of redshift z but have the 
asymptotic values when redshift z becomes large, so at 
the early stage, the VCC models have some properties of 
the quiessence model. The different fate indicated by the 
two models are, if redshift z tends to be —1, the EEoS 
in H power tends to be —1, so the Universe enter a de- 
Sitter phase in the future. However, the EEoS of a power 
law model is always a negative constant, indicating the 
Universe will enter the "Big Rip" phase in the future. 



The cosmological constant problem is still one of the 
serious problem that puzzles the physical world and we 
are still far away to go to its nature. Thus, we expect 
that a more sophisticated combined analysis of various 
observations will be capable of determining the indices 
value of VCC models and revealing more properties of 
the VCC dark energy models. 
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APPENDIX A: SOME RESULTS OF INTEGRALS 



We define a dimensionless function 
E{z) = H(z)/H , 

and integral 

, dz 

I[z) 



(Al) 



(A2) 



la E(z'Y 

For the a power law model, using Eq. (|3ip . we have 

h{z) = r [(l + z)-f+ 1 

(2 - m)Di 

xF 1 (a u f3 ini ,5 1 (l + z) 3 - m ) 

-F l (ai.ft, 71, h)], (A3) 



where 



2 - to 1 

2(3 — m) 2 

8 - 3to . C 3 

71 = 7T75 Z7T' <5l = -77"' 



2(3 -to)' 



C 2 



(A4) 
(A5) 
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and Fx (a, /3, 7, x) represents the hypergeometric func- 
tion. 

For the H power law model, using Eq. (|66p , we obtain 



h{z) = H^(l + z) 



r[(l + z)-^ +1 

(2-m)D| 

xF 1 (ai,/3 1 , 7l ,<5 1 (l + *) 3 - m ) 

-Fi(ai,/?i,7i,^i)]. 



(A6) 



APPENDIX B: DATA ANALYSIS FOR THE 
NUMERICAL FITTING 

We utilize several data sets to constrain the parameters 
of the two power law model. The free parameters in 
these two models are power law index m (or n), current 
value of fractional energy density of dark matter f2 m o, 
and the current value ofh (h = Hg / '100/ km- s -Mpc ). 
However, since the current value of Hubble parameter h 
is always around 0.70 which is determined by current 
supernovae constraints so we marginalize it and plot the 
contour maps of power index m and n versus fl rn o- Our 
data sets include 307 high quality "Union" SN la data, 
baryon acoustic oscillation measurement from the Sloan 
Digital Sky Survey, 9 observational H(z) data and the 
shift parameter from WMAP three years results. 



1. Selected high quality SN la data set 

The first standard candle we use is the type la super- 
novae (SNe la), which is published by Supernova Cos- 
mology Project (SCP) team recently [37]]. This data set 
contains 307 selected SNe la that includes several current 
widely used SNe la data, such as Hubble Space Telescope 
(HST) @,|3|, SuperNova Legacy Survey (SNLS) [MEl 
and the Equation of State: SupErNovae trace Cosmic Ex- 
pansion (ESSENCE) H3. The likelihood function can be 
determined by \ 2 statistics, for the type la supernovae 



2 _ (Pth (parameters; z t ) 



E 



/^cxp) 



r *2 



(Bl) 



where 



\ith (parameters; z) = 5 log <1l{z) + 25, (B2) 



where dh is the luminosity distance which is determined 
by Eq. (TM| 



d L (z) = (1 + z) 



dz 



= H^(l+z)I(z). (B3) 



Baryon Acoustic Oscillation measurement from 
SDSS 



In the large-scale clustering of galaxies, the baryon 
acoustic oscillation signatures could be seen as a stan- 



dard ruler providing the other important way to con- 
strain the expansion history of the Universe. We use the 
measurement of the B AO peak from a spectroscopic sam- 
ple of 46,748 luminous red gala xies (LRGs) observations 
of SDSS to test cosmology [42j, which gives the value of 
A = 0.469(W0.98)-°- 35 ± 0.017 at z B ao = 0.35 where 
n s = 0.95 [43]. The expression of A can be written as 



.4 



1 



(H(z BAO )/H 

V"mO 



ZBAO 



dz 



H(z')/H > 



3 



E(z BAO )z ^BAO 



and the Xbao ^ s 



Xbao 



•A-0.469(n s /0.98) 
. 0.017 



-0.35 x 2 



(B4) 



(B5) 



3. Observational H(z) Data (OHD) 

By using the differential ages of passively evolving 
galaxies determined from the Gemini Deep Deep Survey 
(GDDS) and archival data [44|, Simon et al. determined 
H(z) in the range < z < 1.8 [45j].The 9 observational 
H(z) pieces of data could be obtained from [HI, |4(| and 
they have been have been used to constrain the dark en- 
ergy potential and equation of state [4(| . The \ 2 statis- 
tics for these H{z) data is 
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xohd = 

i=l 



(H (parameters; Zi) — Hi) 2 



(B6) 



4. CMB Data from WMAP three years results 

The CMB shift parameter may provide an effective way 
to constrain the parameters of dark energy models since 
it has the very large redshift distribution and be able to 
constrain the evolution of dark energy very well. The 
shift parameter R which is derived from the CMB data 
takes the form as 



r z 

R = V"mO / 
JO 



dz' 



H{z')/H 



Vernal {zCMb) 



(B7) 



The WMAP3 data gives R 
have 



1.70 ± 0.03 [471, thus we 



XCMB 



/R- 1.70\2 
V 0.03 J 



(B8) 



To break the degeneracy and explore the power and dif- 
ferences of the constraints for these data sets, we use 
them in several combinations to perform our fitting: SN 
+ BAO, SN + BAO + OHD, and SN + BAO + BAO + 
CMB. 
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